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Abstract 

Sensor networks are particularly applicable to the tracking of objects in motion. For such applications, it may not necessary 
that the whole region be covered by sensors as long as the uncovered region is not too large. This notion has been formalized by 
Balasubramanian et.al. [1] as the problem of A^-weak coverage. This model of coverage provides guarantees about the regions in 
which the objects may move undetected. In this paper, we analyse the theoretical aspects of the problem and provide guarantees 
about the lifetime achievable. We introduce a number of practical algorithms and analyse their significance. The main contribution 
is a novel linear programming based algorithm which provides near-optimal lifetime. Through extensive experimentation, we 
analyse the performance of these algorithms based on several parameters defined. 

I. Introduction 

Wireless sensor networks are being increasingly used to due to its vast spectrum of applications including surveillance, security, 
healthcare, asset tracking and environmental monitoring. Such a network consists of large numbers of nodes which are placed 
in a region in order to gather and transmit information about activities in that region. Each node is a low power device 
consisting of sensors and communication units with limited processing power. In general, since sensor networks deployment 
and replacement of nodes in a sensor network is an expensive process, it is of great interest to conserve energy resources in 
the network and extend the lifetime of the network. 

The coverage problem in sensor networks (see [4], [5]) has been extensively studied in literature. The worst-case coverage 
problem may be defined in terms of two related problems: (i) The Maximal Breach problem , discussed in [7] and (ii) The 
Minimal Exposure problem , discussed in [6]. The Maximal Breach problem defines it in terms of the distance of an object to 
be tracked from any sensor and the Minimal Exposure problem defines it by using the sensitivity of a point by a sensor and 
finding the integral of the sensitivities of the points on the path of a moving object to be tracked. In both cases, paths where 
these parameters have lowest values are used to characterize the coverage of a region. 

In many regions, sensor nodes are deployed in large numbers randomly and this generates a lot of redundancy in the network. In 
order to maximize the coverage time, sensors may be turned on and off so that coverage is still maintained. Various techniques 
have been proposed for maximizing the lifetime by allowing certain nodes to sleep (see [3], [9], [10]). In [1], the authors 
investigate the coverage time characteristics of sensor networks and develop algorithms to provide optimal coverage time when 
sensors are allowed to be turned on and off. 

With regard to sparse coverage, most of the work has used the notion of fractional coverage, i.e., only a fraction of the region 
is allowed to be uncovered. The majority of such research involves trying to make the coverage holes occupy as little area as 
possible (see [8], [11]). 

In general, it may be advantageous to model sparse coverage using the diameter of coverage holes instead of area. For 
example, in surveillance applications, it may be necessary to restrict the motion of a target under surveillance. The target 
may be restricted to only move a certain distance before it is detected. In order to solve this problem, the notion of hi-wesik 
coverage was introduced by Balasubramanian et. al. [1]. Recently, it was independently re-introduced by Balister et.al. [2] as 
trap coverage and was investigated for a random deployment of sensors. To the best of our knowledge, [1] was the first work 
to introduce this model of coverage and provided a greedy algorithm for achieving /^-weak coverage. In /^-weak coverage, the 
parameter k, denotes the restriction on the diameter of the holes. 

There are two interesting points in the definition of /^-weak coverage, when compared to full coverage, on top of the much 
longer lifetime. 

1) It is easy to compute, and as our experiments indicate, much easier to give a PTAS (polynomial time approximation 
scheme) to within factor of 1 + e for the coverage holes. 

2) It is easier to compute in a distributed manner. 



Our Contributions: In this paper, we study several practical algorithms for /^-weak coverage, with theoretical guarantees. 
We introduce several algorithms to find hi-wcak coverage including a novel algorithmic correlation between network flow 
algorithms and hz-weak coverage. We provide a Linear Programming algorithm based on this correlation. 

In all settings, it is important to distinguish between the preemptive and non-preemptive cases. In the former, we can turn 
sensors on and off (activate and deactivate) as long as their total active time is at most their original battery life and in the 
latter, once a sensor is activated, it remains active till its battery is depleted. In this paper, we assume that their original battery 
is, for convenience, uniform, but in general is does not have to be the case. We will also normalize such that if a sensor is 
working for time S than its battery loses S%. We also assume that the range of communication is twice the range of sensing. 
So two sensors can communicate if and only if their sensing regions overlap. 

The rest of the paper is organized as follows. Section |Il| deals with the formalization of the problem and the objectives of 



this paper. Section III goes over the techniques employed and the lifetime guarantees achievable. Section IV talks about the 



algorithms to find /^-weak covers together with their scheduling in detail. Section |V] discusses the simulation results and in 



Section VI we provide the conclusion. 

II. Problem Formulation 

The network model consists of a set of sensors S which are placed in a region of interest R known as the Sensor Field . The 
sensing range of each sensor i located at (x^y), is a disk of radius 1 centered at its location and is denoted by Di. Initially, 
every sensor has the same battery life and this is assumed to be 1 unit of time. A point p is said to be covered if p e Di for 
some sensor i. From now on, when we say a sensor i covers a region, it implies that the region is contained in Di. We restrict 
ourselves to square environments of dimensions s x s. 

Definition 1. Let n > 1 be a fixed parameter. We say that a pair of points (p, q) is a n-pair if the distance between them is 
A curve ji is a n-curve is its endpoints form a n-pair A set of sensors C ^ S is a hi- weak cover if every n-curve that is 
contained in R intersects |^ Ds. 

sec 

The problem under consideration is to obtain a collection C = {Ci, C2, . . . , C/c} of subsets of S together with a schedule 

k 

A = {61^62^ ... ^Sk} where each Ci is a hc-wcak cover and C = ^^^i is maximized. The parameter C is termed as the 

i=l 

Lifetime of the system. 

We are also interested in finding a PTAS for the above problem. In this case, we are interested in finding C and A such that 
each Q G C is a ((1 + e)K,)-wQ3k cover. 

In addition, we define the following types of schedules. 

Definition 2. A schedule A is said to be (i) Non-Preemptive if for any Si e A, Si = 1, i.e., each cover Ci is kept ON till the 
batteries are used up, (ii) Uniform if for any Si^Sj G A, Si = Sj, i.e., all covers are kept ON for the same amount of time 
and (Hi) Non-Uniform if for some Si., Sj G A, Si ^ Sj, i.e., all covers need not be kept ON for the same amount of time. 

III. Overview and Lifetime Guarantees 

We first give an upper bound for the optimal lifetime of the system. Then, we proceed to outline the motivation for the 
Grid-Based Techniques and prove some useful properties of these techniques which provides a guarantee of the lifetime of the 
system under these techniques. 

Definition 3. The following definitions are with respect to a set of sensors S. The depth of a point p, denoted by dp is defined 
as the number of sensors in S which cover p. The depth of a region R denoted by dji is defined to be the maximum depth of 
any point p in R. 

Let denote the maximum lifetime of the system under optimal scheduling. The following lemma gives an upper bound for 
Theorem 1. £* = 0(/^ • dn). 

Proof: Let C* = {Q , C|, . . . , C^. } be the covers and A* = {^^ , . . . , } be the schedule which provide the optimal 

k* 

lifetime L*, so that ^ S* — L*. Pick an arbitrary line / and assume that / is the interval between the two most remote sensors 

i=l 



along In each time interval S*, let 5^,85 G C* be two sensors along / such that there is no sensor which covers the portion 
of the line between Sa and 55. The distance between Sa and 55 is atmost k,, otherwise C* is not /^-weak. 

Hence, the minimum number of sensors in C* which intersect / is £ and the total battery life used up for (5* is ■ ^. Now, 
the the total battery at the start of the schedule for the set of sensors intersecting / is atmost f ■ dR. Therefore, 

i=l 

k* 
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L* = 0{k • dR) (1) 

■ 

In this paper, we formulate two techniques in order to find C and A: (i) Grid-Based and (ii) Random Seeds. 
A. Grid-Based 

Consider a tiling T of in which each component of the tiling has diameter of at most hz. Now, it is clear that any subset 
of sensors C that completely covers the boundaries of T is a /^-weak cover for R. The rest of this section investigates the 
lifetime of the system by generating C from such tilings. 

Definition 4. The Total Edge Length of a tiling T = {Ti, T2, . . . , T^}, denoted by TELj- is defined as the sum of the lengths 
of the boundaries of the tiles in T. More formally, if It denotes the length of the boundary of T, then, 



TELr = I 



T 

Ter 



Intuitively, for any tiling T, if ^ 1 then in any cover C for T, \C\ is directly proportional to TELq-. 

We investigate the properties of two tilings: (i) the square tiling and (ii) the hexagonal tiling. From now, whenever we refer 
to a grid, it implies that we are referring to the edges of a tiling. 

Definition 5. A Square Grid for a region R of dimensions L x L is defined to be a set of vertical and horizontal line segments 
{/i, • • • 5 ^m} such that each li is of length L and each square of the implied tiling is of dimensions ^ ^ (^^^ Figure 
^^for examples of such grids). 



Definition 6. A Hexagonal Grid for a region R of dimensions L x L is defined to be a set of line segments {/i, ^2, • • • , ^m} 
such that each li is an edge of a hexagonal tiling superimposed on R where each hexagon is of diameter n (see Figure \2b\for 
examples of such grids). Hence, the length of each U is |. 

Lemma 1. The TEL of the square tiling is 4>/2^. 

Proof: The size of each square is Hence, the number of squares inside an L x L region is Since the perimeter 
of each square is 4^, we have 

since each edge is counted twice. Hence, TEL = 4^/2 — . 



Lemma 2. The TEL of the hexagonal tiling is 

Proof: The side of each hexagon is | since its diameter is k, and its area is • k,'^. Hence, the number of hexagons 
inside an L x L region is -4= • Since the perimeter of each hexagon is 3/^, we have 
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T E L H QxagonalGrid — ' ' 



since each edge is counted twice. Hence, TELnexagonaiGrid = ^l^- ■ 

As input to the algorithm, we are given a collection of such grids Q = {Gi, G2, • • • , G^} and as output, we find a set of covers 
C such that each subset G^ G C is a cover for one of the grids Gj G G- The algorithm iterates over all the grids repeatedly 
until no cover is found for all the grids in Q. The steps of the algorithm are outlined in Figure [T] In essence, there are two 
parts to this algorithm: (i) Obtaining the collection Q and (ii) Finding a cover for a grid G. We describe (i) here and (ii) is 
described in Section |IV] Further, we also describe how to find A in Section [rv] 

Remark 1. We may infer that the two criteria for obtaining maximum lifetime from a collection of grids Q are (i) \Q\ and (ii) 
The number of sensors common between any two grids Gi , Gj ^ Q. If the second criteria is minimized, each cover may be 
activated for longer periods of time implying that the lifetime increases. This can be achieved by selecting grids with lower 
TEL. 



Grid-Based Algorithm 

Generate collection of grids Q based on shifting mechanism, 
repeat 

Find Cover C for grid next G e Q. 

Schedule G, i.e., reduce battery life of sensors in G according to policy, 
if G is tha last grid in Q then 
Restart from the first grid in Q 
until no covers found for all grids in Q. 

Fig. 1. Outline of Grid-Based algorithms. 



1) Finding Q: We describe a method of obtaining Q by considering all possible placements of a tiling in R such that the 
number of sensors which may cover any two tilings Gi and Gj is minimized. Let us examine the square and hexagonal grids 
separately. 

Square Grid: In the case of the square grid, we start with an underlying grid of cells gx g. Let this grid by Gu- Now, consider 
some square grid G such thats its vertices are points in Gu, and some square P in the tiling implied by G. Let be a vertex 
of P. Now, we may obtain other square grids by shifting vp into points of Gu inside the original P along the diagonal of P. 



The grids induced by vp will be the required Q (see Figure 2a). 



Hexagonal Grid: In the case of the hexagonal grid, we start with an underlying triangular grid Gu where each triangle has 
side g. Let each grid G be represented by some hexagon P and let one of the vertices of P be vp. Now, we may obtain other 



hexagons by shifting vp to points of the triangular grid inside the original hexagon P (see Figure 2b). 



In both cases, the value g, which may be termed as Granularity is key to ensuring that the covers are as disjoint as possible. 
In the case of the square grid, if the value g is set to 2, then it can be seen that a sensor can be part of at most two grids Gi 
and Gj. This is because, for any two vertical(resp. horizontal) lines k and Ij part of grids Gi and Gj respectively, no sensor 
can cover both li and Ij . Hence, a sensor can cover two grids Gi and Gj only at an intersection point of lines vertical line k 
and horizontal line Ij from grids Gi and Gj respectively. Similar values can be computed for the hexagonal grid also. 

2) Lifetime Guarantees: Setting the values of ^ to 2 and x respectively, we obtain the following lemmas. 

Lemma 3. If a set of sensors S provides full coverage for R, then the lifetime of the system when using the square tiling to 
find individual covers is at least — when the schedule is uniform. 

Proof: As seen in Section |III-A1| setting the value of ^ to 1 implies that no sensor can be part of two vertical lines in 
two different grids. Hence, each sensor can be part of at most 2 covers, one for which it intersects a vertical line and the other 
when it intersects the horizontal line. The number of shifts when ^ = 1 is given by how many points are on the diagonal of a 
square of the tiling which are at least \/(2) apart. Hence, the number of shifts is Since each sensor can be part of at 

most 2 covers, we may activate each cover for \ units of time. Hence, the lifetime of the system is at least • \ = ■ 

The lifetime of the hexagonal grid can be similarly shown to be ^^^^ . This shows that the guaranteed lifetime of the 
hexagonal grid obtains better lifetime than the square grid. 



3) Approximation Scheme: Let the cover for a particular grid G be G. Consider a line segment / G G and consider the sensors 
Gi (Z G which cover I. Let STRIPi be a strip of width e containing / and equidistant from / on either side and let the 
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(a) Possible Square Grids are shown in Green, Blue and Red 



(b) Possible Hexagonal Grids are shown in Green, Blue 
and Red 



Fig. 2. Shifting Mechanism for Grid-Based algorithms. In both the square and hexagonal cases, an underlying grid is used (square and triangular respectively). 



The points of shifting Vp. 



Vp, c 



are shown inside the polygon P shown with a dashed boundary. The value g is the granularity. 




(a) Approximate Cover for I when I is an internal (b) Approximate Cover for I when I intersects the 

segment. boundary of R. The leftmost sensor is not restricted 

to cover pi but only to intersect the boundary of R. 



Fig. 3. e-approximate cover shown for a line segment I from pi to p2. All sensors in the cover( shown as grey disks) are only restricted to lie within a strip 
of width e which contains /. The requirements are (i) The sensors should be connected and (ii) If / is an internal segment, the points pi and p2 should be 
covered and if / intersects the boundary of R, some sensor should intersect the boundary of R within the strip. 



boundary of RhtVR. In the approximation scheme, Ci need not cover the whole of /. Instead,it is only required to cover 
STRIPi (see Figure |3]). Specifically, the requirements for the cover Ci differ based on whether / intersects the Tr, 

1) Ci should be connected 



2) If / is an internal segment, then Ci must cover both endpoints of / (see Figure 3a). Otherwise, if one or both ends of / 
intersect F^^, then Ci must cover some portion of F^^ within STRIPi (see Figure [3b|). 

This scheme yields covers where there may be curves whose endpoints are at most (1 + e)/^ away from each other. In other 
words, the cover is a (1 + e)/^-weak cover. 

In Section |V| we see that this yields great dividends in lifetime. 
B. Random Seeds 

The random seeds technique involves selecting at random a subset S e S sls seeds and turning off sensors surrounding these 
seeds to form holes as large as possible while maintaning the /^-weak property. The algorithms using which the sensors are 



turned off depends on whether GPS is available and is described in Section IV 



IV. Finding a Cover C and an activation time 5 
In this section, we describe several heuristics to find a cover C for a given grid G and an activation time S for C. 



Definition 7. A Minimal n-weak cover is defined as a n-weak cover C for R such that C\s for any s G C is not a n-weak 
cover 



In the following algorithms, we are interested in finding minimal covers always. This is because, for any cover which is not 
minimal, the redundant sensors would worsen the second criteria described in Remark [T] 

A. Finding a cover for the Grid-Based Algorithm 

In each of the following algorithms, we are given as input a grid G = {/i, ^2, • • • , ^/e}. Let the graph A be the unit disk graph 
of sensors S, i.e., we have a vertex corresponding to each sensor and an edge for every two sensors whose ranges intersect. 
A is unweighted and undirected. 

1) Min-Max Heuristic: The Min-Max Heuristic is described by Balasubramaniam et.al. [1] for full coverage. The idea is that 
at each stage, the sensors are arranged in decreasing order of their battery life. Then, the algorithm checks whether the removal 
of sensor s ^ S (selected in order) violates the cover. If yes, it is made active and otherwise, made inactive. For more details, 
refer to [1]. 

In our case, instead of checking whether removal of a sensor violates full coverage, we check if removal of a sensor violates 
the coverage of the grid G. Let s be the sensor under consideration and let / be a line which intersects s whose endpoints 
are pi and p2- Also, let Ai = (Si^Ei) be the subgraph of A induced by sensors Si which intersect I. Now, removal of sensor 
s may violate coverage only for such lines. Hence, when we remove s, we check if all lines / are covered. This is done by 
checking if the remaining sensors in Si are connected and also whether pi and p2 are covered. In the case of the approximation 



scheme, the second requirement is modified as described in Section III- A3 



Breadth-First-Search Algorithm 



Given a grid G = {/i, ^2, • • • , ^fe} and graph A from S and a paramater e such that < e < 1. 
Let Tji be the boundary of R. 

Let S = 62, • • • , ^m} be the set of possible battery lives of sensors S. 

Perform binary search on B to find the maximum lifetime bj such that a cover is found for G with battery life bj. 
At each step of the binary search when b is the current battery, do the following steps, 
for each line k e G with endpoints pi and p2 do 

Let STRIPi- be the strip of width e • n around li equidistant on either side. 

Let 5'^ C 5* be the set of sensors s such that U STRIPi^ 7^ (j) and B{s) > b. 

Let Ss (resp. Sd) be the set of sensors coverin pi (resp. P2). For the approximation scheme, if pi (similarly P2) 
lies on F^^, let Sg (similarly Sd) be the set of sensors which intersect the portion of F^^ inside STRIPi.. 
Augment Ai. with two vertices s and d and add edges from s to Ss and d to Sd- 
Find the shortest path tt from s to d. 
Ci = {s : s e pi}. 
C=[jCi. 
leG 



Fig. 4. Steps of the Breadth-First-Search Algorithm 
2) Breadth- First-Search (BFS): We describe the BFS algorithm for the approximation scheme with some parameter e such 



that < e < 1. In the exact case, simply set e to 0. F^^ and STRIP are defined as in Section III- A3 Now, for some line 
/ G G with endpoints pi and p2, consider some set of sensors Si which intersect STRIPi. Let Ai be the subgraph of A 
induced by Si. Now, we define two sets of sensors Ss and Sd as follows. If / is an internal segment, then Ss (resp. Sd) be 
the set of sensors covering pi (resp. P2). If Pi (similarly P2) lie on F^^, then, let Ss (similarly Sd) be the set of sensors which 
intersect the portion of F^^ inside STRIPi. Figure [5] shows how Ss and Sd are generated. 

Let B = {^1, ^2, • • • 5 ^m} be the set of possible battery lives of S. We perform binary search on B to find the cover with 
maximum battery life. At an intermediate stage of this search, when b is the current battery, the set Si is found by finding 
sensors intersecting / which have battery life at least b. This is done for each / G G. Then, for each Ai, the augmentation with 
two vertices s and d is performed as described above. Now, all we need to do is find if there is a path from s to d for every 
/ G G (see Figure |5]). If such paths exist, the sensor constituting the paths will form the minimal cover G. See Figure [4] for 
the steps of the algorithm. 

The time complexity of the BFS algorithm is 0(m • ^ • n), where m is the number of unique battery values. This is because 
there are at most m times we may perform the search for the shortest path and at most - lines and each search takes 0(n) 
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Fig. 5. Example of augmented graph and the result of the Breadth-First-Search Algorithm. Here, pi lies on Tfi and p2 is an internal point. Two additional 
vertices s and d are added, s is connected to all sensors intersecting Fji within STRIPi. d is connected to all sensors within STRIPi covering p2- The 
result of the BFS algorithm is the path tt shown as BLUE fat lines. All sensors in tt form the cover Ci. 
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(a) Augmentation of Aq with additional vertices s, d and ji 
and edges from and to fi, from s, to d and between boundary 
sensors intersecting adjacent lines. Also shown are the four types 
of regions in R. 




Bis, 




^k.out 



(b) Replacement of a single sensor Si into Si^in and si^out- sensors 
{si,52,...,Sfc} adjacent to Si have been replaced and the corresponding 
edges are split into two edges from {si^out, S2,out, • • • , Sk,out} to Si^in and 
vice versa. The capacities of these edges is unconstrained. The capacity of the 
internal edge is the battery B(si). 



Fig. 6. Formation of directed graph A^p from Aq. 



time. This is better than the Min-Max heuristic because, for the Min-Max heuristic, the check when each sensor is removed 
may take 0(n) time leading the overall complexity to be O(n^). In practice, the values ^ and m are constant and hence the 
time complexity of the BFS algorithm is only 0(n). 

3) Maximum Flow: Here, we describe a maximum flow algorithm for the square grid case. Given a grid, the maximum flow 
algorithm finds a number of covers C and activation times A. The algorithm has two parts: (i) First, we find the maximum 
flow and (ii) Second, we find individual paths carrying flow where each one is thought of as a cover and the time for the cover 
is the flow through this path. 

Finding Maximum Flow via Linear Programming: Let the sensors which intersect the grid G be Sq and let Aq be the subgraph 
of A induced by Sq- We form a directed graph Alp from Aq as follows. Each sensor Si is replaced by two vertices Si^in 
and Si^out with an edge {si^in, Si^out)- For each edge {si.Sj) of Aq, we have two edges (s^^ont, 5^-,m), {sj^out,Si,in) in Alp. 

The region R is divided into 4 types (See Figure): (i) Vertical Only (marked red in Figure), (ii) Horizontal Only (marked 
green), (iii) Mixed (marked yellow) and (iv) Non-Grid (marked Gray). We define two types of flow in Alp, hence, it is a 



multi-commodity flow problem. On each edge {si^outi Sj,in) ^lp^ we can have horizontal (resp. vertical) flow only if both 
Si and Sj intersect the horizontal(resp. vertical) region. 

Let the horizontal and vertical lines be ordered as Gh = {^i ? ^2 ? • • • ? ^fc} = {^1 7 ^2? • • • 5 ^fc}- Let 5'/^ (resp. Sf^j^) 

denote the set of sensors which intersect STRIPih and the right (resp. left) boundary of and let SJb (resp. SJt) denote 
the set of sensors which intersect STRIPi^ and the bottom (resp. top) boundary of Vr. For each pair of sensors s^, S5 such 
that Sa G S^^ and G S^_^-^ ^ for some i, we add the edge {sa^out-, ^h.in)' Similarly, this is done for sensors in S^_^-^ ^ to 
S^_^2 L ^1^^ similar edges for the vertical lines. Finally, we add three vertices 5, d and ji io Alp and edges 

are added as follows: 

1) For every sensor Sa ^ ^il^ the edge (s, 5a,^n)• 
2) For every sensor G s)^rj., we add the edge {s^^out^d). 

3) For every sensor Sc G S^j^, we add the edge {sc^outilA- 

4) For every sensor Sc G p, we add the edge (/i, Sd,in)- 

See Figure |6] for an illustration of these additional edges. There are no constraints on the capacity of edges between two 
sensors. On each horizontal edge and on edges from s, only horizontal flow is permitted and on each vertical edge and on 
edges to d, only vertical flow is permitted. The internal edge {si^inj Si^out) has capacity equal to the battery life of Si, B{si) 
and the flow permitted depends on which region Si intersects. The purpose of the vertex /i is to convert horizontal flow into 
vertical flow. Hence, edges into /i have horizontal flow and edges outward from /i have vertical flow. For any vertex other than 
s and d, the incoming flow is equal to the outgoing flow and the flow out of s is equal to the flow into d. 

Let the flow on a horizontal (resp. vertical) edge (u^v) be Xu,v (resp. yu,v)- Then, linear program for finding maximum flow 
is given in Figure [7] 

Maximize \^Xss. 



subject to: Wsi ^Xs^^Si 




(2) 




= ^'^ys^,Sj 
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Si 
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^ys,4 
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Si Si 

Fig. 7. Linear Program for finding Maximum Flow 



Finding Invidual Covers: The solution to the linear program gives a flow on each edge m Alp and the maximum flow 
corresponds to a set of covers Cq for G and an activation time 5c for each cover C G Cq- Each cover C corresponds to a 
path from s to d with flow 5c- Hence, C may be found by finding the augmenting path from s io d m Alp. This is done by 
greedily exploring from s along edges containing maximum flow. Once a path pi is found, the corresponding cover C is the set 
of sensors pi. tt is then deaugmented, i.e., the flow along every edge in tt is reduced by the flow along pi. Note that these paths 
may contain loops, in which case, every loop is shortened to get the flow assignment to each edge. Once the deaugmentation 
is complete, the cover may be activated based on the Non-Uniform activation policy described in Section |IV-C| The Uniform 
activation policy does not apply to this method since the optimal solution to the LP gives the activation time for each cover 
in order to get maximum flow. The Non-Preemptive activation policy also does not apply since the optimal solution to the LP 
may cause a sensor to be activated for two different covers in Cc- This process is repeated until no more paths are found. 



B. Finding a cover for the Random Seeds Algorithm 

The random seeds algorithm is designed to operate when GPS is not available, i.e., the sensors do not know their exact location. 
During the course of the algorithm, nodes are labeled as "Boundary", "Deactivated". Initially, all the nodes are unlabeled. The 
algorithm works by randomly picking a set of k seeds and performing Breadth First Search from each of them. The search 



proceeds through unlabeled nodes and stops when a node which is either marked "Boundary" or which is at exactly hops is 
reached, i.e., the search is not continued from such nodes. Also when a node is reached by searches from two or more seeds 
simultaneously, it is marked as "Boundary" and the search is not continued further. At the end, those sensors which are exactly 
hops are marked as "Boundary" and others found during the search are marked "Deactivated". This is then repeated for all 
the seeds. 

In our preliminary simulations, the random seeds algorithm does not achieve comparable lifetime compared to the other methods. 
However, this is the only algorithm which does not require GPS and hence, is worth looking into to make improvements. 

C. Finding an activation time S for a cover C 

We describe three methods to assign activation time 5 for a cover C. These three methods correspond to the three types of 
schedules described in Section |II1 

Uniform: In the Uniform activation policy, each cover is assigned the same amount of time for which to be active. The method 
in which this is assigned follows the description of the Min-Max heuristic in [1]. 

Definition 8. The Load of a sensor, denoted by load{s) is defined as the number of covers in which s participates. 

We assign each sensor s a maximum value for load{s) at the start and at the discovery of each cover, we reduce the battery 
of the sensor by ^. Hence, after m covers, the battery will be depleted. In this manner, we may find the maximum number 
of covers such that the maximum load of every sensor is as assigned. Note that this policy will not be applicable in the case 
of the maximum flow algorithm since the maximum flow algorithm gives a solution in which the durations of the individual 
covers may not be equal. 

Non-Uniform: In this policy, after finding each cover C, we find the sensor s e C which has minimum battery b. Then, we 
reduce the batteries of all the sensors in C by a value d • min B{s), where d is a user-defined parameter termed as the decay 

coefficient . 

Non-Preemptive: In the Non-Preemptive case, once a cover is found, it is activated for the entire battery life of the sensors. 



V. Experimental Results 



In our simulations, sensors are placed in a 50 x 50 region R. R is then divided into cells of size 1x1 and sensors are placed 
randomly inside the cells according to a spatial poisson process of some Intensity . We generated 50 such sets of sensors of 
intensities 1 and 2 and the mean of the lifetime values was calculated. The algorithms Min-Max and BFS algorithms were 
simulated with all three activation policies in order to compare them with the square grid. In the case of the hexagonal grid, 
we simulated the BFS algorithm with uniform activation policy to compare its performance with the square grid. The LP was 
simulated with non-uniform activation policy for the square grid. Each of these simulations was performed with = {10, 20} 
and e = {0, 0.1, 0.2, 0.3}. The granularity of the shifts was set to g-\-e, where g is set to the values from Section III-A1[ since 
we want to capture as disjoint sensors as possible with each shift. Figure [8] shows four example covers for the square grid for 
e = {0, 0.1, 0.2, 0.3}. It can be seen that as e increases, the sensors deviate further and further from the lines of the grid. 

In Figures [9j [T0| [TT] and 12 Min-Max , BFS and LP denote the corresponding algorithms executed on the square grid and 
HexGrid denotes the BFS algorithm on the Hexagonal Grid. Figures [9] and 10 show the lifetime vs e for = 10 and k, = 20 
respectively when the intensity of the poisson process is 1 and Figures [TT] and 12 show the same for intensity 2. In each of 
the figures, it may seen that the lifetime values increase sharply first followed by a decrease of slope in the curve. The net 
lifetime is, however, improved greatly over the lifetime when e = 0, i.e., when the algorithms are exact. This decrease in slope 
may be attributed to the fact that, as e increases from ei to 62, the ratio of the number of sensors considered for ei to those 
considered for 62 gets lower. This, in turn, seems to contribute directly to the decrease in slope. 

Figures [9a| [9a| and [9c] and similarly, in Figures [9] [TO] [TT] and [12] show the distribution of lifetime vs e in the case of the three 
activation policies. In all cases, the BFS seems to perform comparable to the Min-Max algorithm with better performance for 
the Non-Preemptive case. Hence, we are able to achieve similar performance with lower communication complexity. It may 
be seen in the Non-Uniform case (Figures [9b] |lQb[ [TTb] and |12b| ) that the LP algorithm performs much better than the others. 
This is in line with expectations since the LP provides close to optimal scheduling for covers found for a single grid, i.e.. 
Finally, in Figures [9a| 1 1 Qa[ [TTa| and 1 1 1 a| we see that the hexagonal grid offers a very big jump in lifetime over the square grid 



algorithms, again, in line with the expectations considering the guarantees provided in Section III-A2 




(c) e = 0.2 (d) e = 0.3 

Fig. 8. 



VI. Conclusion and Discussion 

In this paper, we have investigated the problem of /^-weak coverage in sensor networks. The idea of using tihngs to find 
A^-weak covers was introduced with hfetime guarantees provided for the Square and Hexagonal tilings. A novel LP algorithm 
was introduced to find covers for such tilings and its performance was compared against two other algorithms and was shown 
to achieve significantly better lifetime. The hexagonal grid was compared against the square and was also shown to achieve 
significantly better lifetime. Finally, an algorithm for finding covers in case GPS was not available was introduced. 
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(c) Non-Preemptive 
Fig. 9. Lifetime vs e when k, = 10 and Intensity = 1. 



(c) Non-Preemptive 
Fig. 10. Lifetime vs e when k = 20 and Intensity = 1. 




(c) Non-Preemptive 
Fig. 11. Lifetime vs e when k, = 10 and Intensity = 2. 



(c) Non-Preemptive 
Fig. 12. Lifetime vs e when k, = 20 and Intensity = 2. 



